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Abstract 



By the classical Martingale Representation Theorem, replication of random vectors can 
be achieved via stochastic integrals or solutions of stochastic differential equations. We in- 
troduce a new approach to replication of random vectors via adapted differentiable processes 
generated by a controlled ordinary differential equation. We found that the solution of this 
replication problem exists and is not unique. This leads to a new optimal control problem: 

find a replicating process that is minimal in an integral norm. We found an explicit solution 
00 ' 

^T) ' of this problem. Possible applications to portfolio selection problems and to bond pricing 



models are suggested. 
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o5 '■ 1 Introduction 

By the classical Martingale Representation Theorem, random variables generated by a Wiener 
process can be represented via stochastic integrals. This means that it is possible to find a 
Ito process such that the terminal value matches a given random vector at a fixed terminal 
time. This result leads to the theory of backward stochastic differential equations and has many 
applications in Mathematical Finance. 

We introduce a new approach to replication of random vectors via adapted differentiable 
processes generated by a controlled ordinary differential equation. We found that the solution of 
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this replication problem exists and is not unique. Therefore, an optimal control problem arises: 
to find a replicating control process that is minimal in a certain norm. We found an explicit 
solution of this problem in linear quadratic setting. Possible applications to portfolio selection 
problems and to bond pricing models are suggested. 



2 The control problem 

Consider a standard probability space (0, T, P). Let w(t) be a standard d-dimensional Wiener 
process which generates the filtration Tt = cr{w(r) : < r < t} augmented by all the P-null sets 
in T; we assume that w(0) = 0. 

Let / € L 2 (Q, Tt, P; R- n ) be a random vector. 

Let a £ R ra , A 6 j^x™ ^ e given, and let b E R nxn be a non-degenerate matrix. Let 

r : [0, +oo) — > R, nxn be a measurable matrix valued function. 

Consider the following linear quadratic optimal stochastic control problem 

fT 

Minimize E / u(t) T T(t)u(t) dt over u(-) (2.1) 
Jo 

subject to *W-^) + K<), te(0,T) ^ 
x(0) = a, x(T) = f a.s. 

This is a stochastic control problem with equality type constraints on the value of the plant 
process that holds almost surely. This problem is related to control problems for backward 
stochastic differential equations (BSDEs) with fixed terminal value for the plant process. In 
BSDEs control setting, a non-zero diffusion coefficient is presented in the evolution equation for 
the plant process as an auxiliary control process. The first problem of this kind was introduced 
in [7]. Our setting is different: a non-zero diffusion coefficient is not allowed. Similar problems 
were introduced in [8], [9]. 

2.1 Admissible T and u 

For p > 1 and q > 1, we denote by Lp* m the class of random processes v(t) adapted to Tt with 
values in R,™ xm such that E I J Q \v(t)\ q dt\ < +oo. We denote by | • | the Euclidean norm for 
vectors and the Frobenius norm for matrices. 

Let g : [0, T) — > R be a given measurable function such that there exist c > and a € (0.5, 1) 
such that 

< g(t) < c(T - t) a , (2.3) 
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git)- 1 <c(l + (T-t)- a ), 



t€[0,T). 



An example of such a function is g(t) = 1 for t < T — Ti, <?(i) = (T — t) a for t > T — Ti, where 
Ti 6 (0, T] can be any number. 

Let {/ be the set of all processes from L^ 1 such that 



We consider U as the set of admissible controls. 

We assume that T(t) is a measurable matrix valued function in R raxn , such that T(t) = 
g(t)G(t), where G(t) > is a symmetric positively defined matrix such that the matrices G(t) 
and G(t)" 1 are both bounded. By the definitions, E u(t) T T(t)u(t)dt < +oo for u EU. 

Restrictions on the choice of T(t) = g(t)G(t) mean that the penalty for the large size of u(t) 
vanishes as t — > T. Thus, we do not exclude fast growing u(t) as t — > T such that u(t) is not 
square integrable. 

2.2 Optimal control u 

By the Martingale Representation Theorem, there exists a unique kf G LV^ g such that 




(2.4) 



By the definition of * , it follows that, for u € U, 




(2.5) 




(See, e.g., Theorem 4.2.4 in [14], p. 67). 

We assume that there exists r E (0, T) such that 



esssupE|fcj (£)| 2 < +oo. 



Let 






By Lemma 1 from [9], it follows that k^(-) € 
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Theorem 2.1 Let 

-p(t) = 22(0) -1 (E/ - e AT a) + I k^s)dw(s) 

Jo 

and 

u(t) = T(t)- 1 b T e AT ^m- 

Then this u belongs to U and it is a unique optimal solution of problem i2.1\ )- {KlQ) in the class 
ueU. 

In [8], a related problem was considered for a simpler case when it was required to ensure 
that x(T) = B{f\T e } for some 9 <T. 

3 Applications to finance 

Replication on the basis of Martingale Representation Theorem is the main tool in modern 
Mathematical Finance. The presented above approach to replication on the basis of Theorem 
12.11 can also be applied to problems arising in finance. Some possible applications are listed 
below,including optimal cash accumulation policy and modeling of the bond prices. 

3.1 Optimal cash accumulation policy 

Consider a model where a there is a risky asset with the price S(t) which is a random continuous 
time process with positive values adapted to a Wiener process w(t). We assume that 

dS{t) = S{t)[a{t)dt + a(t)dw{t)}, 

where a(t) and <r(t) are some J^-adapted bounded processes such that a(t) > C a.e., where 
C > is a constant. 

Let terminal time T > be given. 

Assume that an investor wishes to accumulate gradually an amount of cash that allows to 
purchase a share of this risky asset at time T. Let u(t) be the process describing the density 
of the cash deposits/withdravals at time t G (0, T), such that u(t)At is the amount of cash 
deposited during the time interval (t, t + At), for a small At > 0. Assume that it is preferable 
that the cash flow will be as smooth as possible. 

Let us assume first that the bank interest rate is zero, for both loans and savings. In this 
case, the total amount of cash at the terminal time will be u(t)dt. 
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Theorem 12.11 can be applied now for n = 1, A = 0, b = 1, / = S(T). In this case, 

J s 

If a(t) =0 then Theorem 12.11 ensures that the process 

u{t) = T(t)- 1 \R(0)- 1 S{0) + [ Ris^dSit)] (3.1) 
L Jo i 

is such that 

f u(t)dt = f a.e. (3.2) 
Jo 

Moreover, this u(t) is optimal in U in the sense of the optimality criterion from Theorem 12.11 
i.e., E T(t)u(t) 2 dt is minimal. 

Let us consider a more general model where the bank interest rate is r > 0, for both loans 
and savings. In this case, the total amount of cash at the terminal time will be e r ^ T ~^u{t)dt. 

Theorem 12. II can be applied now for n = 1, A = r, b = 1, / = S(T). In this case, 



R( 8 )= [ Q(t)dt, Q(t) = e 2r ^T(t)- 1 . 

J s 



If a(t) = then Theorem 12.11 ensures that the corresponding process (|3.ip is such that 

rT 

e r{T -^u{t)dt = f a.e. (3.3) 

Again, this u(t) is optimal in U in the sense of the optimality criterion from Theorem 12. 1\ i.e., 
E Jq T(t)u(t) 2 dt is minimal. 

If a(-) ^ 0, then conditions (j3.2|) are still satisfied for u(t) defined by ()3. 1 j) but the value 
E Jq T(t)u(t) 2 dt is not minimal over u anymore. Instead, Eq T(t)u(t) 2 dt is minimal, where 
Eq is the expectation defined by an equivalent probability measure Q such that S(t) is a mar- 
tingale; we will call it a martingale measure. This still means that deviations of u are minimal 
but in a different metric. It can be also noted that the definition of the class U for the original 
measure has to be adjusted for the new measure Q, with the expectations E replaced by Eq. 

Let us consider a modification of the cash accumulating problem where the accumulated 
cash amount has to be a given proportion of the excess achieved by the equity at the terminal 
time. This problem arises for a writer of a naked or a partially naked call option. To cover this 
case, it suffices to apply Theorem 12 . 1 1 with / = cmax(5(T) — K, 0), where c > is the prescribed 
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proportion. We assume that a(t) is non-random and that the bank interest rate is r > 0. In 
this case, it is well known that 

t-T 



f = H(S(0),0) + J -^(S(t),t)dS(t). 



Therefore, k f (t) = %Z(S(t),t)a(t)S(t) and Eg/ = if(5(0),0). By TheoremEU is ensured 
for this / with 

u(t) = cT(t)- 1 R(0)- 1 H(S(0),0)+ / i2(s) _1 -3-(5(t),*)rf'S(*) » 

where H(x,t) = Eg {max(5*(T) — K, Q)\S(t) =x}. This expectation is under the martingale 
measure Q again; if a(t) is constant, it can be calculated by the standard Black-Scholes formula. 
In addition, Eg T(t)u(t) 2 dt is minimal among all u such that (I3.2p is satisfied. 

Another possible choices of / in this setting may include / = cmin <e [ ,T] S(t) (a given pro- 
portion of the minimum achieved by the equity during the time period [0, T]), or / = y Jq S(t)dt 
(a given proportion of the average equity value over the time period [0, T]). Here K > and 
c > are some constants. In these cases, u(t) also can be represented explicitly for non-random 
a(t). 

The model described above can also be applied to the problems of optimal dividend flow 
selection. In particular, it can be applied to the setting where the manager of a firm with the 
capitalization S(t) wishes to pay dividends during the time period [0, T] such that the total 
payoff Jq u(t)dt over this interval will be, say 5% of the equity S(T) at time T. The typical 
approach is a barrier criterion of dividend payments or analysis of ruin times; the methods are 
usually based on dynamic programming (see, e.g., [5], [11]) and the bibliography here). Theorem 
12.11 leads to a new approach to this problem. 

3.2 Modelling of the bond prices 

Consider continuous time bond pricing model for zero coupon bonds. Let B(t,T) be the bond 
price at time t for the zero coupon bond with payoff $1 at time T, where T > t. Let r(t) be the 
short rate. We assume that the process r(t) is J^-adapted. Here Tt is the same as above; it is 
the filtration generated by a Wiener process. 

We assume that the probability measure P is a measure used for the pricing such that, for 
a given process r(t), 

r-T 



Let B(t, T) be the discounted bond price denned as 

B(t,T) = exp(- J r(s)ds^B(t,T). 

In particular, 

(■T 



B(T, T) = exp^— J r(s)ds 



Usually, a model evolution of r(t) is being suggested first and then the distributions of B(t,T) 
and B(t,T) are derived. Theorem 12.11 gives an alternative: to model first the distribution of the 
random variable £ = B(T, T), and derive the evolution low of r(t) from the distribution of £. It 
appears that it can be done for a quite wide class of random variables £ with the values in (0,1) 
such that / = — log£ satisfies conditions of Theorem 12.11 Since £ € (0, 1), we have that / > 
a.s. By Theorem 12.11 there exists an adapted process r(t) such that 

f = J r(s)ds, £ = exp^— J r(s)dsj. 

Moreover, the process r(t) can be selected to be optimal meaning that it has minimal deviations 
(in the sense of the optimality criterion from Theorem 12. lj) . 

This approach can be extended on the case of a bond market where there are bonds with 
different non-random maturity times T/., k = 1,...,N, > T^+i for all k. Assume that we 
are given J-j> h -measurable random variables with values in (0,1). Let fi = — log£i, /2 = 
log£i — log£2r--> fk = log£fc — log£fc_i. Applying Theorem 12.11 modified for the positive initial 
times, we obtain that there exists an adapted process r(t) such that 

fi= r(s)ds, f k = / r(s)ds, k = 2,...,N, 



JO 

and 

(■T k 



£ k = exp(^-J r(s)dsj, k = l,...,N. 



This leads to a bond market model such that 

B(T k ,T k ) = ^k, k = l,...,N, 

for an arbitrarily chosen set such that the corresponding random variables f k has final 

second moments and that condition (|2.6p is satisfied for T = T k and / = f k . As we had 
mentioned already, the conventional approach is to select a model for the process r(t) first and 
then to derive B{T k ,T k ). This possibility to start with models for B(T k) T k ) is established here. 
This could give new opportunities for modeling of bond prices. 
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4 Proof of Theorem 12.11 



For the sake of completeness, we give below the proof of Theorem 12. 1| this proof follows the 
proof of Theorem 1 from [9] . 

Clearly, equation (|2.2p gives that 

x (t)= f e A{t - s ^bu{s)ds + e At a. (4.1) 
J o 

By ([53]) . x(T) G L 2 (n,J r T ,P;R n ) for any u G U. 

Let the function L(u,fi) : U x L 2 (Q, Tt, P; R n ) — > R be defined as 

L(u, n) = ~E f u{t) T T(t)u(t) dt + E^ T (/ - x(T)). 

For a given //, consider the following problem: 

Minimize L(u, fi) over u EU. (4-2) 

This problem does not have constraints on terminal value x(T). Therefore, it can be solved 
by usual stochastic control methods for the forward plant equations. We solve problem (|4.2p 
using the so-called stochastic maximum principle that gives a necessary condition of optimality; 
see, e.g., [6]- [T] , [E]-[T3], [H]-[l6]). For our problem (|4.2]h all versions of the stochastic 

maximum principle from the cited papers are equivalent and can be formulated as the following: 
if u = € U is optimal then 

m T bu^(t) - ^(t) T T(t)u^t) > m T bv - ^v T T(t)v (4.3) 

for a.e. t for all v € R" a.s., where ip(i) is a process from L r 2 l ^ such that 

dij){t) = -A T ip(t)dt + x(t)dw(t), 
4>(T) = M, 

for some process % S -^2 2™- (See, e.g., Theorem 1.5 from [3], p. 609). The only solution of the 
backward equation for ip is 

= e^V-Q^t), n{t) = E{n\T t }. (4.4) 

Necessary conditions of optimality (|4.3p are satisfied for a unique up to equivalency process 
u = defined as 

u„(t) =T(t)- 1 b T iP(t)- (4-5) 
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Let us show that £ U for any \i. We have that 

e ^ \ Ufl (t)\dt\ < CiE ^ T irct)- 1 !!^*)!^ 



<C 2 sup E|^(i)| 2 / 0(i) -1 dt < +00. 
te[o,T] io 



In addition, 

E [ T g(t)\u^t)\ 2 dt<C 3 B /" s(t)|r(t)-V(t)| 2 di <C 4 E ! g{t)- l \^{t)\ 2 dt 
Jo Jo Jo 

r T 

< C 4 sup E|//(t)| 2 / g(ty l dt < +00. 

t6[0,T] Jo 

Here > are constants defined by A, 6, n, and T. Hence € £/. 

Clearly, the function L(u, /u) is strictly concave in u, and this minimization problem has a 
unique solution. Therefore, this u = is the unique solution of (|4.2p . 

Further, we consider the following problem: 

Maximize L(u^,^i) over \i G L 2 (fi, Jr, P; R n ). (4.6) 
For u = u^, equation (|4.ip gives 

X (T) = [ T e^-'hu^dt + e AT a. 







Hence 



L(n M , /x) = -E / ^(t) T r(t> M (t) tft - E^" 1 " / e^^bu^dt - Bfi T e AT a + E^ T /. 
2 io Jo 



We have that 

rT 



E/i T / e^^bu^dt = E^ T / e A ( T - i )6r(t)- 1 6 T ^(t) ( it 
io io 

rp rp 

= E / Li T / e A{T - t) bT(t)- 1 b T e AT{T - t) fi(t)dt = Efj7 [ Q(t)n(t)dt 
Jo Jo 

= E / Li T Q{t)^(t)dt = E [ E{fi T Q{t)^(t)\T t }dt = B [ n(t) T Q(t)n(t)dt. 
Jo Jo Jo 

The fifth equality here holds by Fubini's Theorem. Further, we have that 

e/ u fl {t) T r(t)u f ,(t)dt = B [ (r(t)- 1 6 T ^(t)) T r(t)r(t)- 1 6 T V'(t)^ 

io io 

= E f tP(t) T bT{t)- 1 b T ^j{t) dt = E / (e AT(T ~*V(i)) T fcr(t)- 1 6 T e AT ( T - t V(t) dt 
io io 

= E / ^(t) T e A ( T -*)6r(t)- 1 6 T e AT(T -*V(i) dt = E [ fi{t) J Q(t)n(t) dt. 
Jo Jo 



It follows that 

1 r T 

L( W ) = Bfi T (f - e AT a) - -E J fi(t) J Q(t)/i(t) dt. 
By the Martingale Representation Theorem, there exists G such that 

p = fl + k,j,(t)dw(t), 
Jo 

where p, = E/i. It follows that 

E^ fi(t) T Q(t)fi(t)dt = B (ji + J kp(8)dw(a)J Q(t)(fL + k IJL {s)dw{s) S j dt 
= p T Q{t)pdt + B J (J k^dw^s)) Q(f) (J kp(a)dw(8)J dt 

= / \/ Q ^ dt )f i + J V N k^s)dw(syj Q{t) (J k^{a)dw{a)\ dt 
= /Z T i?(0)/2 + I dtE I k fl (s) T Q(t)k^s)ds = p T R(0)p + E [ ds [ k^s) 1 Q^k^s) dt 

JO Jo Jo Js 

= pJ R(0)fi + E / ^(s) T ii(s)A; M (s)ds. 

JO 

We have used Fubini's Theorem again to change the order of integration. Similarly, 

B fJ T f = p T f + n f kp(t) T k f (t)dt, 

Jo 

and E// T e AT a = p T e AT a. It follows that 

L(« MJ p) = p T (f - e AT a) - lp T R(0)p - ^E F ^(r) T R(r)k^T) dr + E F k^(t) T k f (t) dt. 

1 1 Jo Jo 

Clearly, the maximum of this quadratic form is achieved for 

p> = R^THJ -e AT a), k^t) = RQty^flt). 
This means that the optimal solution p of problem (|4.6p is 

p = R(0)- 1 (f -e AT a)+ f kJt)dw(t). 







Let u(t) and be defined by (I4.4j) - (l4.5p for fj, = ju, i.e., 2 = up. By Lemma 1 from [9], it 
follows that kfi(-) € ££,2 d - It follows that 



Jo 



E / \i,„U)\ 2 dt < +00. 
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It follows that sup te r T ] E|/I(t)| 2 < +00. 

We found that sup^ inf u L(u, fx) is achieved for (u, /2). We have that L(u, fi) is strictly convex 
in u € U and affine in fi € £2(^5 J 7 , P, R n ). In addition, L(u, fi) is continuous in u £ -^i? 1 gi ven 
// G ^2(^1 -T 7 ! R n )i and L(u, /u) is continuous in ^ E L2(fi, J 7 , P, R n ) given u £ U. By 
Proposition 2.3 from |10| . Chapter VI, p. 175, it follows that 

inf supL(ti,^) = sup inf L(u,fi). (4-7) 

Therefore, (u,fi) is the unique saddle point for (|4.7p . 

Let Uf be the set of all u(-) 6 £/ such that (12. 2 p holds. It is easy to see that 

rT 



I f 

inf -E / u(t) T T(t)u(t)dt = inf sup L(u, //), 
ue(7/ 2 Jo «et/ ^ 



and any solution //) of (14 .7p is such that u £ Uf. It follows that u £ Uf and it is the optimal 
solution for problem (|2.ip - (|2.2|) . Then the proof of Theorem 12.11 follows. □ 
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